Spectra and binding energy predictions of chiral interactions for ^Li 
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Using the no-core shell model approach, we report on the first results for based on the next-to- 
next-to-leading order chiral nuclear interaction. Both, two-nucleon and three-nucleon interactions 
are taken into account. We show that the p-shell nuclei are sensitive to the subleading parts of the 
chiral interactions including three-nucleon forces. Though chiral interactions are soft, we do not 
observe overbinding for this p-shell nucleus and find a realistic description for the binding energy, 
excitation spectrum and radius. 
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I. INTRODUCTION 

We are currently experiencing rapid progress in our 
understanding of nuclear properties. This is triggered by 
two major developments. 

On one hand, we have increasingly powerful supercom- 
puters, for which new and very efficient algorithms have 
been developed to solve the nuclear many-body problem. 
We are now able to solve the Schrodinger equation for 
realistic nuclear interactions for p-shell nuclei including 
also three-nucleon forces (3NF's) 0, This is a ma- 
jor advance in itself, because it is becoming more clear 
that reliable predictions for many nuclear observables, 
binding energy and spectra can be obtained from phe- 
nomenological nuclear two- (NN) and three-nucleon (3N) 
interaction models. Especially, for the NN system, these 
have reached a high degree of sophistication and describe 
the NN data up to pion production threshold perfectly 
0, 0, 1^, For the binding energies of p-shell nuclei, 
the structure of the 3NF's turned out to be significant, 
leading to improved models of these forces engineered 
to describe a wide range of light nuclei accurately 0- 
This tool is of increasing importance to determine, e.g. 
reaction rates for astrophysical processes @, which are 
experimentally not accessible or properties of nuclei with 
large neutron excess [J. 

On the other hand, there was a great deal of progress 
in our understanding of how chiral perturbation the- 
ory (ChPT) can be extended from purely pionic or the 
nucleon-pion system (for a review see e.g, . 13) to sys- 
tems with more than one nucleon |lOl IllL Il2j |. In this 
approach, one makes use of the explicit and spontaneous 
breaking of chiral symmetry to systematically expand the 
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strong interaction in terms of a generic small momentum. 
Thereby, the NN interactions, the SNF's and also ttN 
scattering are related to each other. The chiral symme- 
try and the pattern of its breaking are not systematically 
taken into account by today's phenomenological interac- 
tions, except that in all of them the longest range part of 
the potential is generally the one pion exchange interac- 
tion. Therefore, though QCD is believed to be the theory 
of the strong interaction for the energies of relevance in 
nuclear systems, we are not able to perform confirming 
tests using the traditional forces. This will be possible 
using nuclear forces based on ChPT. It will be especially 
important to look at subleading parts of the interaction, 
which include the 27r exchange NN and 3N forces. Here, 
many of the relations between NN, 3N and ttN interac- 
tions become apparent. Therefore, finding signatures of 
the 3NF's is an important aspect of current research on 
this issue. 

In the past, the effects of 3NF's have been studied using 
the phenomenological models. In nucleon-deuteron (Nd) 
scattering above 60-100 MeV lab energy, it was found 
that predictions for some (polarization) observables de- 
pend on the 3NF model used, but not on the NN interac- 
tion chosen jlj]. These observables are excellent candi- 
dates to pin down the structure of 3NF's and, therefore, 
outstanding laboratories to study the important sublead- 
ing parts of the chiral interactions. Consequently, a lot 
of experiments were triggered, which provide an impor- 
tant set of data to probe these models 0, 0, 0, 0| . 
First comparisons for chiral interactions to the data were 
performed and, in general, agreement was found in the 
energ y ra nge, where chiral interactions are expected to 
work [131 • 

These data, however, are manifestly isospin T = 1/2 
and, therefore, not suitable to probe the 3NF in the 
isospin T — 2>/2 channel. But it was argued that the 
T = 3/2 component is very important to describe p-shell 
nuclei, which mainly motivated the new 3NF terms in the 
Illinois series of 3N interactions [7f. This clearly shows 
that the spectra and binding energies of light nuclei are 
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complementary to nd scattering and provide indispens- 
able information on the structure of 3NF's. It is there- 
fore of utmost importance to predict properties of p-shell 
nuclei based on the chiral interactions including these 
3NF's. 

Chiral interactions are low momentum interactions. 
The unknown short distance part of the force is absorbed 
in a tower of contact terms. Necessarily, we need to reg- 
ularize the interactions with cutoff functions. It turns 
out that one obtains a decent description of the NN data 
using rather small cutoffs. This is advantageous, because 
the resulting interactions are soft and convergence is gen- 
erally faster than the convergence with traditional mod- 
els. But this has also been criticized in the past. The 
experience with traditional models indicated a need for 
rather hard cores not only in the NN interaction, but 
also in the 3NF to prevent strong overbinding in sys- 
tems beyond the s-shell. This is a surprising result, be- 
cause it means that we cannot properly separate high 
energy degrees of freedom from the low energy part. Be- 
cause the high momentum tail of traditional interactions 
is strongly model dependent, a sensitivity to this tail will 
induce a strong model dependence of the results. Indeed, 
the binding energies of nuclei do depend on the chosen 
NN interaction model. It is believed that this model 
dependence is strongly reduced, once 3NF's consistent 
with the NN force are added (in fact it was proven that 
any NN interaction can be augmented by a 3NF so that 
the description of the 3N data is identical to that of a 
second phase-equivalent NN force , thereby removing 
the model dependence). Similarly, one can expect that a 
consistent combination of NN and 3N force will lead to 
a reasonable description of light nuclei even if the high 
momentum tail is completely missing. An important as- 
pect of this work will be to confirm this expectation by 
an explicit calculation. 

Because of the strong non-locality of today's chiral in- 
teractions, to the best of our knowledge, there is only one 
reliable many-body technique to solve the Schrodinger 
equation for p-shell nuclei for these interactions includ- 
ing the 3NF's. This is the ab-initio no-core shell model 
(NCSM) d 13, ISH. In this paper, we use the NCSM 
to predict spectra and binding energies for ^Li. These 
results will be presented along with calculations for the 
s-shell states for '^H and ^He, which were used to deter- 
mine unknown parameters of the 3NF. Thep-shell results 
are then predictions. 

The manuscript is organized as follows. In Section ITll 
we start with a brief summary of the status on the re- 
search on chiral interaction models, mainly to explain 
the force model entering in the calculations. The NCSM 
approach is introduced in Section lllll to keep this paper 
self-contained. Then, in Section IIVI we discuss our nu- 
merical results. Conclusions and an outlook will be given 
m Section El 



II. CHIRAL INTERACTIONS 

Microscopic nuclear structure is generally based on NN 
interactions, which describe the NN scattering data up 
to TT-production threshold perfectly. Usually, such in- 
teractions are based on the l7r-exchange, which is aug- 
mented by a more or less phenomenological short range 
part. It turns out that binding energies based on these in- 
teractions are model dependent and too small P, . 
3NF's are clearly necessary. 

However, defining consistent combinations of NN and 
3N forces is not a trivial task, and rarely are both based 
on a single underlying theory. In practice, the NN inter- 
action models are usually augmented by a 2tt exchange 
3NF in or extended versions thereof 0- Then 
parameters are mostly adjusted to reproduce the bind- 
ing energies of s-shell nuclei or nuclear matter density. 
Though this can be very successful in the description of 
nuclear properties , it lacks a solid theoretical founda- 
tion. The relation to QCD is completely lost. Processes 
of other strongly interacting systems, like the ttN system, 
are not quantitatively related to these forces, though the 
basic mechanisms involve these particles. 

Since lattice simulations for nuclei are not realistic for 
the foreseeable future, effective field theory is an appeal- 
ing theoretical foundation on which to build. In its sim- 
plest form, an effective field theory was formulated, that 
explicitly takes only nucleons into account. It is accom- 
panied by a power counting in terms of powers of r/a, 
where r is the effective range of the interaction and a the 
large scattering length of the NN system (for a recent re- 
view see pdj*). This property is realized in very different 
physical systems, e.g. ^He clusters and atoms close to a 
Feshbach resonance. It allows one to identify universal 
properties. Also, it is possible to identify nuclear proper- 
ties, which are dominated by the large scattering length 
of the NN interaction. An example is the correlation 
between '^H and "^He binding energies, known as Tjon- 
line, which is naturally explained using this approach 
|27| . While these results are interesting in themselves, 
they are by construction not that interesting for the re- 
lation of QCD and nuclear physics. The approach will 
only describe observables, which do not reflect any spe- 
cific property of QCD except the large scattering length. 

To sense the leading role of QCD, another effective 
field theory is the tool of choice: ChPT. The new phys- 
ical input is the approximate chiral symmetry of QCD, 
which is known to be spontaneously broken. Based on 
this an effective field theory can be formulated that in- 
volves nucleons and tt's explicitly. The Goldstone-boson 
character of the tt's (related to the spontaneous breaking 
of chiral symmetry) guarantees that amplitudes in purely 
pionic and the ttN system can be expanded in powers 
of a small, typical momentum Q. The explicit break- 
ing of chiral symmetry due to the finite quark masses is 
expanded similarly in terms of the tt mass l^^. Elec- 
tromagnetic effects on the quark level and quark mass 
differences have also been taken into account by an ex- 
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pansion of isospin symmetry breaking contributions in 
terms of the fine structure constant and the quark mass 
differences |29ll3nj. 

Obviously, the NN system is different. The deuteron 
bound state clearly shows that a perturbative treatment 
is impossible. Weinberg reahzed that this is understood 
by an enhancement of reducible diagrams in two- or 
more-nucleon systems and proposed to expand the ir- 
reducible diagrams perturbatively in terms of small ex- 
ternal momenta, but then sum all these diagrams to in- 
finite order u sing a Lippmann-Schwinger or Schrodinger 
equation [sj- In other words, Weinberg suggested 
to obtain nuclear potentials from ChPT. This naturally 
explains that NN interactions are much more important 
than 3NF's and that these are more relevant than higher- 
body interactions [33|. The approach was quantitatively 
confirmed [33. IssH S q j. For the NN system, terms up to 
order momentum Q"^ have been taken into account nowa- 
days ^2Iie and lead to a description of NN data, which 
is comparable to the one of the phenomenological inter- 
actions. The approach to sum all irreducible diagrams 
to infinite order is usually called "Weinberg counting". 
Other schemes to deal with the non-perturbativity have 
been proposed and are formally more consistent, be- 
cause the renormalization is understood analytically. But 
the formal problems of the "Weinberg" scheme, which 
were pointed out in [s^, can be circumvented by small 
cutoffs of the order of 500 MeV for the regularization of 
the Lippmann-Schwinger equation. Under this condition, 
"Weinberg power-counting" is numerically more success- 
ful than the scheme of ;39i as was shown in [Z^l . Indirect 
confirmation of the Weinberg approach is the favorable 
agreement of parameters extracted from ttN scattering 
[m and an NN phase shift analysis li^ . Therefore, 
for the rest of this paper, we will rely on the Idaho-N3LO 
interaction .37j based on "Weinberg counting" for the NN 
interaction. It takes terms up to order Q'* into account 
and is able to describe the NN scattering data below tt- 
production threshold perfectly. We stress that this can 
only be a first step, because a complete analysis of the 
chiral interactions requires further order-by-order analy- 
sis to confirm convergence. 

Very soon the leading 3NF's, which appear at order 
were derived They consist of a 2tt exchange 

term, which is v ery similar to the widely used Tucson- 
Melbourne 3NF z5\ and, additionally, there are contact 
interactions with and without Itt exchange. In the orig- 
inal paper, five such contact terms were believed to be 
independent. Fortunately, it turned out later, that only 
two of them remain independent once the identity of the 
nucleons is taken into account. In this simplified form 
the 3NF's were applied to Nd scattering, the '^H and ''He 
bound states 0. For completeness, the explicit expres- 
sions for these forces are given in Appendix IXI 

For the 27r-exchange term, the power of this system- 
atic scheme to derive consistent NN and 3N interactions 
comes into play. All vertices also occur in NN interac- 
tion diagrams. This completely determines the strength 
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Fettes et al. [45j -1.23 -5.94 3.47 

Biittiker et al. [41j -0.81 -4.70 3.40 

Rentmeester et al. [42| -0.76 -4.78 3.96 

Entem et al. [37j -0.81 -3.20 5.40 

Table I: Various determinations of the strength constants a 
of the 3NF. All values are in GeV"^ 



and form of this part of the 3NF. As outlined in the ap- 
pendix, there are three less well known strength constants 
in this part of the 3NF: ci, C3 and C4. They determine the 
contribution of various spin-isospin structures to the off- 
shell ttN scattering amplitude entering the 3NF. Table HI 
summarizes the results of recent determinations of these 
constants. The most fundamental determinations were 
done in Ref. |4ll li^j ] using a fit to ttN scattering data. 
For this system, ordinary ChPT without resummation is 
used. The possibility to extract these constants from ttN 
data makes the strong link of the nuclear interaction and 
other strong interaction processes in ChPT explicit. The 
rather scarce set of data for the ttN system induces rather 
larg e uncertainties. This has partly been circumvented in 
using dispersion relations. Therefore, the determina- 
tion of Ref. |4l| is usually considered to be more accurate. 
Therefore, it is nice to see that this result compares very 
well with the determination of where the q's were 
extracted from a partial wave analysis of NN data. For 
the Idaho-N3LO interaction, the authors of ^3 preferred 
their own determination for C3 and C4 based on a fit to 
phase shifts in high partial waves, also shown in the ta- 
ble. It slightly deviates from the two previous values and 
lead to an intense debate on the correct way to extract 
these constants from NN data. We cannot contribute to 
this discussion, but will use the q's of simply be- 
cause we want consistency of the NN and 3N forces. The 
exact magnitude of these constants should be of less im- 
portance, because they only appear in subleading terms 
of minimal order Q^. However, this issue has to be kept 
in mind, when we compare our results to the experimen- 
tal values, especially for observables, where 3NF's are 
important. Please note that we will combine an NN in- 
teraction of order Q'* to a 3NF of order Q'^ only. Because 
we will not study the convergence of our predictions with 
respect to the chiral expansion in this work, we think that 
this combination is justified. It is the state-of-the-art for 
chiral interaction models. 

For the remaining two contact 3NF's, two more 
strength constants enter: cu and ce (see the appendix 
for the definition). In Ref. 18], it was shown that the 
determination of these constants is possible by a fit to the 
■^H binding energy and the Nd doublet scattering length. 
This is only one of the possibilities to determine these 
constants. For the Idaho-N3LO force, we therefore de- 
cided to fit to the ^H and ^He binding energies. For our 
study here, this has two advantages. Firstly, a very accu- 
rate fit is possible, because both binding energies are well 
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Figure 1: Correlation of the two dimensionless strength con- 
stants Co and ce of the 3NF. All combinations of cd and ce 
describe the '^H binding energy equally well. 
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Figure 2: Eigenvalues rj of the Yakubovsky equations for *He 
depending on cd- ce is chosen according to Fig. 77 = 1 
corresponds to describing the ''He binding energy exactly. 





CD 


Ce 


3NF-A 


-1.11 


-0.66 


3NF-B 


8.14 


-2.02 



Table II: Combinations of cd and ce, which describe the H 
and ■'He binding energies. 



known experimentally. Secondly, it makes a comparison 
of our results to calculations with the phcnomenlogical 
model AV18/Urbana-IX more meaningful, because this 
combination describes both binding energies reasonably 
well, but fails for some p-shell spectra. The fit was per- 
formed using Faddeev- Yakubovsky calculations. For a 
detailed description of the technique, we refer to [2^ . 
Using the '^H binding energy, we established a relation 
between cd and ce shown in Fig.^ For all these combi- 
nations of cjj and ce, Idaho N3L0 in combination with 
the 3NF predict a binding energy of 8.482 MeV in 
agreement with experiment. These combinations then 
enter in a calculation of the Yakubovsky eigenvalue 77 
for ^He using the a particle binding energy of 28.3 MeV 
as trial energy. Our theoretical prediction agrees with 
28.3 MeV, if rj is equal to one. From the results shown in 
Fig. 121 we can read off that two combinations oi cd and 
Ce describe and ''He equally well. For later reference, 
we call them "3NF-A" and "3NF-B" . The numerical val- 
ues are listed in Table ITU 

We did not perform complete calculations for the bind- 
ing energy for ^He for each combination of c_d and ce- 
However, generally, one can expect from the ranges of 
•q results that the binding energy would change by ap- 
proximately 2 MeV in the range of cd's considered. This 
sets a scale for the contribution to *He. For com- 
pleteness, we also note that the pp force was augmented 
by the Coulomb interaction for point like protons. No 



further electromagnetic contribution was taken into ac- 
count as they are expected to be of less importance. We 
also do not include the first relativistic corrections to the 
3NF, which some people do count as an order con- 
tribution. We cannot discuss this issue here. We only 
claim that, whatever counting you prefer, the relativistic 
corrections will be smaller than our , because we use 
cutoffs around 500 MeV, which are much smaller than 
the nucleon mass. Then, because no four-nucleon forces 
enter at order Q'^, our interaction is complete up to this 
order. Now we need to extend our calculations to the 
p-shell. For this, we will give a brief introduction to the 
NCSM in the next section. 



III. NCSM APPROACH 

Usually shell-model calculations exploit the possibility 
to assume an inert nuclear core. Taking only the valence 
nucleons as active particles clearly has the advantage 
to keep the number of many-nucleon states manageable, 
but, so far, it also means that effective shell model inter- 
actions have to be used, which cannot be related directly 
to nuclear interactions, as they have been developed for 
few-nucleon systems. The "ab-initio no-core shell model 
(NCSM)" or, more simply, the NCSM is different 
Because all nucleons are taken to be active, in principle, 
the same interactions can be used for the NCSM as for 
traditional few-body calculations. However, for realistic 
nuclear forces, the wave function is hardly described by 
combinations of single particle states. Using harmonic 
oscillator (HO) states, which allow for an easy separa- 
tion of the center of mass (CM) motion from the internal 
degrees of freedom, one encounters the additional prob- 
lem to describe the the exponential long range tail of 
the wave function. Again effective interactions become 
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necessary to achieve convergence, with a small number 
of HO states taken into account. These effective interac- 
tions are systematically related to the "bare" interactions 

mill- 

The scheme described in Ref. |i47J is based on unitary 
transformations of the Hamiltonian |48l | , which decouples 
the model space from the complete Hilbert space describ- 
ing the quantum mechanical system. Such transforma- 
tions have already proven to be useful in a wide variety 
of nuclear physics problems. They allow us to decouple 
parts of the Fock space, such as the two nucleon space, 
from the rest of the Fock space ji^ . Also a decoupling 
of low and high momentum components has been per- 
formed using the same approach |^ . In many-body 
problems, it has not only been used for shell model calcu- 
lations, but, recently, could be exploited for the solution 
of nuclear problems in the hyperspherical harmonics ba- 
sis 

The starting point of all NCSM calculations is a non- 
relativistic ^-body Hamiltonian, which includes NN and, 
important here, 3N interactions 

A ^2 A A 

^^ = E£^+ E ^^^+ E ^'^^ (1) 

i=l i<j=l i<j<k=l 

Pi are the momenta of the A particles with masses m^. 
In the following, we assume that the interaction does not 
change the type of particle and that the total mass Ma — 
J2i T^i is well defined. Adding the CM HO potential, one 
can rewrite the Hamiltonian as [S^l 

i<j — 1 ^ 
A 

+ E ^^j-fc (2) 

i<j<k=l 

with the CM position R = -pj rriiri. 

In the following we would like to establish an unitary 
tranformation of the Hamiltonian, which decouples two 
parts of the Hilbert space. A rather small finite model 
space P and the rest of the Hilbert space Q. The pro- 
jection operators on the two spaces are also called P and 
Q. They fulfill the relation Q = l-P. 

Okubo 48] showed that the unitary operator can be 
derived from an operator oj, which fulfills the following 
relations: 

cj ^ Q uj P (3) 

= Q e-"^ H^l e+'^ P 

= Q{l-w) H^2 (l+uj) P (4) 

The transformation e'^ already decouples the subspaces 
P and Q, but it is not unitary. Unitarity can be achieved 
using the antihermitian operator S = arctanh(aj — w^) 



and defining the unitary operator U = e"*" . Using the 
properties of S and lu, one finds the explicit expression 

U = e^ (l + Louj^ +uj^ Lo)-^/^ (5) 

The problem is reduced to finding the operator lo. For 
that, one defines a set of dp eigenstates |n) of the Hamil- 
ton operator H^, where dp is the dimension of the model 
space P. Defining uj \n) = Q |n), it is easy to show that 
the decoupling condition Eq. Q holds. This completely 
defines the oj operator for sets of eigenstates |n), for which 
the P \n) are linearly independent. In practice, for res- 
onable choices of the eigenstates, this condition is always 
fulfilled. Using Eq. Q the action of ui on \n) reads 

(agin) = ^{agl w \ap){ap\n) 

a p 

= E^"QI ^ l"p)^apn (6) 
a p 

a enumerates the basis states. The indices P and Q 
restrict enumeration to the respective subspace. With 
the help of the inverse M = 7V~^, the matrix elements of 
ui are found [55| easily in a non-iterative scheme 

(agl UJ \ap) = ^(agln) M„ ap (7) 

n 

As elaborated in [s^, there is always as set of states 
|n), for which the inverse of N exists. In practice, one 
chooses the dp lowest lying eigenstates of H^. Because 
(1— cij + cjT)|n) has only P-space components, the same 
holds for U^n). Matrix elements of states \n) are then 
exactly reproduced in the model space 

(n'l O \n) = {n'\U PU^ O UP U^\n) (8) 

Because we are only interested in the P-space matrix 
elements of Og/ / , we can make use of the simplification 

C/P = (1 + cj) (1+cjt t^)-i/2p (9) 

and write 

{ap\O,ff\0p)= ^ {ap\il + cu^ u;)-'/^\a'p) 

a'p P'p 

J2^nc',M„'t3'p {n\ O In') 

n 7i' 

{P'p\{l + u;^ u;r'/^\Pp) (10) 

Here we made use of the relation 

{ap\{l+uj^) O {l+uj)\(3p) 

= E ^'^ncpMn'Pp {n\ O |n') (11) 

n n' 

There are two special cases of interest of Eq. (Illll . The 
first, O = 1, is useful to determine the P space renormal- 
ization operator (1 -I- since 

(«p|(l+C^t) (l+c^)|/3p) 
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{ap\{l+uj^uj)\(3p) 



(12) 



The second, O — H^Jj makes the evaluation of the effec- 
tive Hamiltonian especially simple 

{aplHefflPp) = {ap\il+u;^) H^{l+u;)\f3p) 



(13) 



with En the eigenenergies of the eigenstates \n). These 
equations are exact. They reduce the problem to one in 
a very small model space P. However, because solutions 
of the problem are necessary to obtain the effective op- 
erator, they do not help to solve the eigenvalue problem 
in this form. 

The NCSM concept becomes useful by approximating 
the effective interaction in an a-body cluster approxi- 
mation. There are two obvious constraints on this ap- 
proach. Firstly, to keep the approximation controllable, 
we require that the effective interaction approaches the 
bare interaction for model spaces P — > oo, and, secondly, 
we expect best convergence, when the nuclear part of 
the interaction acts in the a-body cluster with the same 
strength as in a free a-body system. Then, under the as- 
sumption that only b = 2-body bare forces act, one needs 
to solve the cluster problem for 



a 



1=1 \ 



i<j = l 



E (14) 



2Ma 



Note that in this cluster Hamiltonian, the strength of the 
HO two-body potential depends on the original A. This 
induces a confining meanfield potential for the cluster, 
which simply reflects that the cluster is embedded in the 
nuclear environment. It has the nice property that we can 
base our effective interaction on an infinite set of cluster 
bound states. No scattering states need to be considered. 
The effective interaction for 6-body interactions in the 
a-body cluster approximation for the ^-body system is 
defined as 



V, 



eff 



eff 




(15) 



where Ti^f / is derived from an Hamiltonian operator like 
given in Eq. H14|l . The effective interaction enters into the 
A-body problem with a different weight, that insures that 
the bare interaction is recovered for large model spaces 
I2II21I 



^eff 




4=1 

I {A-ay.{a-by. 
{A-b)\ 



E 



V, 



A,a,b 
eff,ii...i 



(16) 



In this form it becomes apparent that we have a double 
counting problem, when both, 6 = 2 and 5 = 3 interac- 
tions enter the cluster Hamiltonian. This can be solved 
as outlined in and requires two effective interactions. 



The first one, V^j'i'^^- , is based on the Hamiltonian in 
Eq. The second one, ^eff'if^ ia' based on 



Hi 



i -t2 




i<j^l 

a 

- E ^^fc 

i<j<k=l 



(17) 



which additionally takes into account the 3NF. The com- 
bination 



Pi , rrii ir 



eff Z^\2mi 2 

1=1 \ 

(A-a)!(a-3)! 



e//,ii...i 



(A -3)! 

[A - a)\{{A - 3)!(a - 2)! -{A- 2)!(a - 3)!) 



{A-i)\{A-2)\ 



eff, 11- 



(18) 



ii<...<j:a = i 



then fulfills the constraints for our effective interaction. 
Specifically for a = 3, we get 



eff - 




rrii ^2 



il<...<ia = l 

3- A ^ 



v„^;:'f^. (19) 



il<i2...<ia = l 



A~2 ^ ^eff,ii...ia 
ii<...<ia = l 

This interaction will be used in the calculations in the 
next section. It will enhance the convergence so that we 
can obtain meaningful binding energies for p-shell nuclei. 
In the following, we also use the a = 3 approximation, 
when only "bare" NN forces are taken into account, i.e., 
when Vijk — 0. 



IV. RESULTS 

The main subject of this section will be a careful analy- 
sis of the numerical accuracy and convergence properties 
of the results. The cluster approximation induces a size- 
able n dependence, which is driven by the mismatch of 
the long range behavior of the cluster states and the wave 
function for the nucleus. The former ones being confined 
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by a HO behave like Gaussians, whereas the latter ones 
decay exponentially for large distances. Because renor- 
malization of the operators is generally weak for the large 
distances this cannot be corrected by the renormal- 
ization of the interaction, but needs to be taken care of 
by an optimization of the for the problem and by a 
rather large model space, which includes enough basis 
states that the tail is also well approximated. Therefore, 
to get reliable results, calculations for different model 
spaces have to be performed for a range of fi's to iden- 
tify the optimal choice ,58]. 

The model space size is characterized by the maxi- 
mal number N of HO excitations included. The lowest 
^Li HO configuration has four nucleons in the n = 
HO states and three in the n = 1 shell. Therefore, 
the minimal HO energy is Afmin = 3. The A-body 
calculations include states up to Nmax ~ ■N'min + N 
|69| . In this work, we will perform calculations up to 
N = Q including configurations up to Mmax = 9. The 
calculations were performed in a basis of single particle 
Slater-deterimants (m-scheme) using the Many-Fermion- 
Dynamics code 59, 60]. Since we employ HO single par- 
ticle states and all Slater determinants up to Nmax are 
included in the model space, the center of mass part can 
be exactly separated at the end. For iV = 6, the P-space 
dimension was 663527 requiring approximately 3-10® ma- 
trix elements for the A-hody Hamiltonian. The effective 
interactions are based on a = 3 clusters. For each model 
space, we include three-body states up to the maximal 
excitation of three particles in the ^-body configurations. 
For our largest calculation, these are three-body states up 
to A3 = 9 and, for the smaller calculations, this number 
decreases correspondingly. This implies that the effec- 
tive interactions are different for each model space size. 
The results for smaller model spaces do not possess upper 
bound character. 

To solve the cluster equation, we use the Lanczos- 
scheme. The usage of a Jacobi basis for this problem 
guarantees the translational invariance of the effective 
interaction. This is also mandatory to exclude unphys- 
ical center of mass excitations from the model space, 
which would also complicate the selection of the proper 
states \n) for the definition of the lo operator. The clus- 
ter equations are solved using the "bare" interactions for 
Q-spaces up to A/3 = 48 for three-body angular momenta 
up to J3 = I and TVs = 40 for J3 = | to ^ . The largest 

dimension for the A3 = 9 P-space is 157 for = | 
and three-body isospin = i making the same amount 
of bound state solutions necessary in this channel. Alto- 
gether, for all channels between J3 = ^ and J3 = we 
required 2352 states. Convergences for the cluster states 
has been checked using their binding energies. We take 
the charge dependence of the NN interaction into account 
by averaging the NN isospin t = 1 interactions according 
to the prescription given in '61*] . This averaging was per- 
formed for total isospin T = 1/2 ^Li states. All states we 
will discuss here are T = 1/2 states. The Coulomb inter- 
action is not included into the calculation of the effective 




Q. [MeV] 

Figure 3: (Color online) Q. dependence of the binding energy 
for the ^Li J'^ = | ground state using the NN interaction 
only. Results for iV = to A'' = 6 model spaces are shown. 



interaction, but is taken into account as an additional 
"bare" force during the solution of the A-body problem. 
This does not change any results, but will facilitate cal- 
culations for T — 3/2 states later, which can now be well 
approximated using the same effective interactions. 

In Figs. 13 ^ and |S1 we show the Vl dependence of 
the ground state binding energy for ^Li with the Idaho- 
N3L0 interaction only and the combinations with 3NF-A 
and 3NF-B, respectively. The binding energy exhibits a 
distinct minimum for all model spaces and interactions 
shown. The position of this minimum shifts to smaller f2 
for larger model spaces. This shift becomes smaller for 
the larger spaces. Also, the position depends on the in- 
teraction model. This behavior has been observed for all 
realistic interactions and makes the extraction of binding 
energy results more difficult than for the simpler Min- 
nesota type of interactions. Ref. 58] investigates this is- 
sue in detail. There, it was shown that taking the binding 
energy of the minimum for different model spaces results 
in a clear pattern of convergence. It was also shown that 
the results obtained in this way are consistent with re- 
sults based on the "bare" interaction. Therefore, we also 
follow this procedure. The inserts in the figures show the 
convergence of the binding energy based on this idea for 
model spaces with A = 0, 2, 4 and 6. Very clearly we 
observe convergence for Idaho-N3LO without 3NF and 
with 3NF-A. For 3NF-B, the changes in the step A = 4 
to A = 6 are somewhat larger. They are, however, small 
enough to justify that the results are accurate to the or- 
der of the few hundred keV's or one MeV. This will be 
sufficient for a discussion of the quality of chiral interac- 
tions later. 

The results for the excitation energies are shown in 
Figs, ini to Fortunately, the O dependence is rather 
mild and the different excitation energies are very similar 
for the A = 2, 4 and 6 model spaces. The vertical lines 
indicate the 17 of the minima for A = 6 in Figs. 13 to 
|31 Our final extraction for the excitation energies is the 
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13 14 15 
n. [MeV] 

Figure 4: (Color online) Q dependence of the binding energy 
for the ^Li J'^ = | ground state using the NN interaction 
combined with 3NF-A. Results for iV = to = 6 model 
spaces are shown. 
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n. [MeV] 

Figure 6: (Color online) fl dependence of the excitation en- 
ergy of the lowest states of ^Li using the NN interaction only. 
Results for TV = to TV = 6 model spaces are shown. 
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n [MeV] 

Figure 5: (Color online) Q dependence of the binding energy 
for the ^Li — ^ ground state using the NN interaction in 
3NF-B. Results for TV = to TV = 6 model spaces are shown. 
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Figure 7: (Color online) Q dependence of the excitation en- 
ergy of the lowest states of ''Li using the NN interaction com- 
bined with 3NF-A. Results for TV = to TV = 6 model spaces 
are shown. 



TV = 6 result at this value. We also checked the 
dependence of the binding energy for the excited l/2~ 
state. Since the position of the minima is very similar 
to the ones for the 3/2~ state, we are confident that the 
excitation energies for the shown 3/2^, 1/2^, 7/2^ and 
5/2 ~ are accurate. 

Our results for the binding energy are summarized 
and compared to other calculations and the experiment 
m Table Uni Besides our new results, we show results 
for phenomenological interactions from GFMC calcula- 
tions. An earher NCSM results for AV8' and TM' gave 
35.8 MeV. It was based only on a single value of hU, i.e. 
14 MeV, and on three-body cluster states up to TV's = 28. 
For a more reliable binding energy prediction, a complete 
study of the fl dependence is certainly required as well as 
an increased number of three-body cluster states. Pre- 
liminary results clearly show this trend. Note that we 
use As = 48 in the present study. Because of this, we 
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Figure 8: (Color online) fl dependence of the excitation en- 
ergy of the lowest states of ^Li using the NN interaction com- 
bined with 3NF-B. Results for A = to A = 6 model spaces 
are shown. 



model 


Egs [MeV] r [fm] 


NN only 


34.6 2.15 


w/ 3NF-A 


38.0 2.11 


w/ 3NF-B 


36.7 2.23 


AV8'+TM' [63, 64] 


40.5(1) 2.2(1) 


AV18+Urbana IX 


37.5 2.33 


AV18+IL2 [7, 65j 


38.9 2.25 


Expt. [66, 67] 


39.2 2.27 



Table III: Comparison of the ground state binding energy re- 
sults Egs and the point proton rms radius r for ^Li (for N = 6 
model space and Q at the minimum of the energy /Q, curve) 
for chiral interactions and several phenomenological combina- 
tions to the experimental value. The experimental rms radius 
is corrected for the finite size of the proton assuming a proton 
charge radius of 0.81 fm~^. 



omit the earlier NCSM result here and compare to the 
GFMC resuhs. 

As expected, ^Li is underbound for the NN interac- 
tion only. Both 3NF's provide more binding. However, 
in both cases the final binding energy result is still short 
of the experiment by 1.2 and 2.5 MeV, respectively. In 
view of the general expectation that strong repulsion at 
short distances is required to avoid a collapse of nuclei, 
we are encouraged by the slight underbinding. Clearly, 
we do not observe any sign of overbinding, though neither 
the NN force nor the 3NF have a strong repulsive core. 
Both interactions are very soft, yet the binding energies 
are reasonable. For phenomenogical, pure 27r exchange 
3NF's, hke the Tucson-Melbourne (TM') 0, overbind- 
ing already sets in for ^Li. So far, it was believed that 
only the addition of a repulsive core, like in the Urbana- 
IX and Illinois models 0, 13 1 '^^^ '^^''^ ^^^^ overbinding 
problem. Here, we show that the additional structures of 
chiral 3NF's also prevent overbinding. We note that the 
description of the binding energy is best for 3NF-A. 

The final results for the excitation energies are sum- 
marized in Fig. O All combinations of the interactions, 
Idaho-N3LO alone, with 3NF-A or 3NF-B, do predict the 
right ordering for these states. The splitting of the 3/2^ 
and 1/2" states is small. The agreement with this experi- 
mental splitting seems to be superior for 3NF-A. Because 
the splitting itself is very small, this might be acciden- 
tal. More significant deviations of the predictions are ob- 
served for the 7/2~ and 5/2~ states. Both, the position 
of this multiplet and the splitting is strongly affected by 
the 3NF's and the agreement with the experimental re- 
sults is clearly best for 3NF-B. This is in contradiction to 
the binding energy, which was best described for parame- 
ter set 3NF-A. For a better clarification of this situation, 
a study of the dependence on the choice of the q's is 
mandatory. This issue is not as simple as it seems at 
first glance, since the same q's enter the NN interaction. 
So, in principle, a refitting of the NN interaction would 
be necessary. We postponed this issue to a later study, 
where this constraint might be loosened to get more in- 
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Figure 9: (Color online) Dependence of the excitation energy 
of the lowest states of ^Li on the interaction. Results with the 
NN interaction only and with the 3NF-A and 3NF-B included 
are compared to the experimental value. 



sight into the Ci dependence. Finally, we note that part of 
the deviations for the binding energy might be attributed 
to our numerically uncertainty. Therefore, the results for 
^Li favor the parameter set 3NF-B, because of its supe- 
rior description of the spectrum, but further calculations 
will be required before we can conclude that this result 
also holds for heavier nuclei. 

As we discussed before, we consider one main result 
of this first calculation to be that the very soft chiral 
interactions do not predict a collapse of nuclei with in- 
creasing number of nucleons. Of course, this also needs to 
be confirmed for the densities. Therefore, we also show 
results for the point proton radii for the 3/2" state in 
Figs. ^1 to El For this quantity, we found an interest- 
ing dependence on H. and N. For all model spaces, the 
radii decrease with increasing f2. This exactly reflects the 
behavior of the basis HO states. Obviously, there is no 
correlation to the binding energy result. It is comforting 
to note that the slope of this dependence decreases with 
increasing TV. Also, we observed the smallest dependence 
on N for the fl values around the minima of the bind- 
ing energy curves, which are again indicated by vertical 
lines. Therefore, we have good indications that the best 
extraction of the radii are obtained for exactly these 
values for the largest possible model spaces. 

Table Unl also shows the radii obtained in this way and 
compares them to the experimental values and other cal- 
culations. We obtain very realistic results for the chiral 
interactions, which are comparable to the phenomeno- 
logical models Urbana-IX, Illinois and TM'. We do not 
find any indication that soft, chiral interactions fail to 
saturate nuclear systems with a realistic binding energy 
and density. We note that for 3NF-B we even observe 
an increase of the radius though it provides additional 
attraction. 
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Figure 10: (Color online) dependence of the point proton 
rms radius of the ground state of ^Li using the NN interaction 
only. Results for TV = to A'' = 6 model spaces are shown. 
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Figure 11: (Color online) f2 dependence of the point proton 
rms radius of the ground state of ^Li using the NN interaction 
combined with 3NF-A. Results for A'' = to A?" = 6 model 
spaces are shown. 
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Figure 12: (Color online) SI dependence of the point proton 
rms radius of the ground state of ^Li using the NN interaction 
combined with 3NF-B. Results for iV = to iV = 6 model 
spaces are shown. 



V. CONCLUSIONS 



We have presented the first microscopic calculation of 
the low lying ^Li states based on a full chiral interaction, 
including NN and 3N forces. Our results are based on the 
NCSM. We showed that predictions are possible with an 
accuracy that allows one to discuss the basic properties 
of the nuclear force. The results confirm that binding en- 
ergies and spectra depend on the structure of the 3NF's. 
This means that the results for different sets of param- 
eters for the 3NF are different though both models de- 
scribe the and ^He nuclei with the same accuracy. In 
this first application, we restricted ourselves to changes 
in the strength of the short range, or short-range long 
range pieces (D- and E-terms) of the 3NF, which are not 
related to the NN interaction at the order Q^. With this 
constraint, we have already obtained reasonable results 
for both, the binding energies and the spectra, but we 
could not describe both simultaneously within the accu- 
racy, which we expect from our many-body method. 

On the other hand, we know that the strength of 27r 
exchange in the NN and 3N interactions is not well de- 
termined yet. Therefore, it is necessary to establish a 
measure for the dependence of the ''Li and other p-shell 
results on this strength. Since, in principle, this can only 
be done by changing the NN and 3N forces, it is not a 
straightforward and simple task and, probably, for first 
estimates, one will need to loosen the strict consistency 
of NN and 3N forces somewhat. This will be the subject 
of a forthcoming study. 

Here, we were mainly interested to establish that soft 
chiral interactions result in realistic binding energies and 
densities. We therefore also investigated the point pro- 
ton radii and described how these can be extracted from 
NCSM wave functions. The results support that soft, low 
momentum interactions predict reasonable binding ener- 
gies and densities. In the future, this conjecture needs 
to be confirmed for heavier systems. In this respect, a 
study of is a high priority. Results employing the 
TM' 3NF and the AV8' NN force indicate an increasing 
over binding going from A = 7 to A = 10 [g^- It needs 
to be confirmed that the chiral interactions can again 
avoid overbinding. At the same time, has proven to 
be very sensitive to the spin-orbit interactions including 
the 3NF's. Usual NN interactions, hke Idaho-N3LO and 
AVI 8 HjQjEII, predict a wrong ordering of the 3+ ground 
and 1+ excited state. The Illinois type of 3NF's correctly 
change this ordering and also the TM' 3NF affects this 
positively, in contrast to the Urbana IX 3NF. Therefore, 
the predictions of the complete chiral interaction with 
the 3NF are especially interesting for this nucleus. But 
also for this system, the dependence on the Ci's of the 27r 
exchange should be studied. 
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Following the notation of |68'|, the 2tt exchange part 
reads 



U. 



(fe):27i 



ijk 



E 



( 9A 



2\2FJ {q, 2 + ml){q, 2 + ml) ^i'^ * ^ 



(Al) 



where qi is the momentum of the pion exchanged between 
nucleons i and k and 



4cim2 2c3 ^ ^ 



(A2) 

The strength of these terms is completely determined by 
the Ci constants discussed in the main text. The two new 
terms read 



Appendix A: EXPLICIT EXPRESSIONS FOR THE 
CHIRAL 3NF'S 

In this appendix, we would like to summarize the ex- 
plicit form of the chiral 3NF. The explicit expressions 
were first derived in The 3NF consists of terms 

with the three different topologies. The first one is the 
usual 27r exchange. The operator structure is identical to 
the well known Tucson-Melbourne interaction ^25]. Then 
there are two new topologies, which contain also contact 
interactions (see e.g. 18j). In the original work, five inde- 
pendent terms with such structures are derived. In [l8l |. 
it was realized that the number of independent terms is 
smaller once the antisymmetry of the 3N states is taken 
into account. Then only one term of the Itt exchange 
type (D-term) and one term of the pure contact type 
(E-term) are independent. 



y-(fe):l7r _ _ yA 



9A CD cTj ■ qj 



and 



(A3) 



8i^2 F^A^ qj^ + ml 



(fc) ;cont 



ijk 



IE 



2 /-^ P4A ^'^ 



{tj ■ Tk) (A4) 



In these terms two new strength constants cd and ce 
appear. The values, which we determine as described 
in the main text from the binding energies of '^H and 
^He, are based on the choice A^ = 700 MeV. For the 
axial vector coupling constant, we use qa — 1.29. -FV = 
92.4 MeV is the weak pion decay constant and nijr ^ 
138.03 MeV the averaged pion mass. As outlined in |l8l |. 
we regularize the 3NF using a cutoff. Throughout this 
paper, we fix the cutoff to A = 500 MeV. 
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